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1 .  INTRODUCTION 


This  study  was  motivated  by  the  existence  of  a 
unique  explosive-driven-implosion  facility  at  the 
Institute.  Stable  focussed  implosions  can  be  pro¬ 
duced  in  this  reliable  facility  yielding  extremely 
high-pressure  and  temperature  conditions  at  the 
focus.  This  facility  has  been  used  successfully 
for  the  generation  of  hypervelocity  projectile 
flight  and  intense  shock  waves  as  well  as  the 
production  of  diamonds  from  graphite  and  neutrons 
and  i  rays  from  deuterium  fusion  [Ref.  1] .  The 
performance  of  the  implosion  chamber  was  investigat¬ 
ed  extensively  [Refs.  2-8],  In  all  these  investi¬ 
gations  it  was  assumed  that  effects  of  viscosity 
and  heat  conductivity  were  negligible.  However, 
it  is  important  for  a  complete  assessment  of  per¬ 
formance  for  the  implosion  chamber  to  evaluate  the 
effects  of  the  viscous  boundary  layer  on  the  im¬ 
plosion  focus.  It  is  necessary  to  consider  these 
aspects  of  flow  in  several  stages.  The  first  stage 
concerns  the  induced  boundary  layer  behind  the 
gaseous  detonation  wave  in  a  hydrogen-oxygen  mix¬ 
ture  and  is  the  subject  of  the  present  report. 

In  general,  all  detonations  or  explosions  fall 
into  two  categories:  structure-independent  and 
structure-dependent .  For  the  former,  the  chemical- 
reaction  zone  is  very  narrow  immediately  behind  the 
wave  front  and  its  structure  has  no  effect  on  the 
characteristics  of  the  flow  field  behind  it.  For 
the  latter  the  relatively  long  chemical-reaction 
zone  and  its  structure  dominate  the  properties 
of  the  flow  field.  Both  types  of  detonation  were 
studied  extensively  [Refs.  9-17].  In  this  study 
of  2H2*02  detonations  only  the  structure- independent 
wave  applies. 

Mirels  and  Hamman  [Ref.  18]  made  the  first 
attempt  to  solve  che  laminar  boundary  layer  behind 
a  strong  blast  wave  moving  in  a  power-law  path. 

They  used  a  series-expansion  method  which  was  also 
employed  later  by  Chen  and  Chung  [Ref.  19].  This 
method  is  only  applicable  in  the  region  near  the 
wave  front.  Sichel  and  David  [Ref.  29]  made  a 
modification  of  the  shock-tube  boundary-layer  solu¬ 
tion  by  Mirels  [Ref.  21]  and  Hartunian  et  al  [Ref. 

22]  to  calculate  the  heat  transfer  to  the  wall 
behind  a  detonation  wave  with  the  assumption  that 
the  pressure,  temperature  and  velocity  behind  the 
detonation  remain  constant  and  are  equal  to  the  C-J 
values.  These  assumptions  are  restrictive  in 
describing  the  actual  physics  of  the  problem. 

Later  on,  numerical  solutions  were  obtained  by  Liu 
and  Mirels  [Ref.  23]  for  different  one-dimensional 
flows  (planar,  cylindrical  and  spherical)  over  the 
entire  flow  region. 

In  previous  analyses  [Refs.  18,  19,  23]  for 
boundary- layer  flows  behind  strong  blast  waves, 
four  major  assumptions  were  made  in  order  to  sim¬ 
plify  the  analysis:  (1)  The  explosion  wave  moves 
with  a  power-law  path  in  the  time-distance  plane. 
Consequently,  similarity  solutions  exist  for  the 
inviscid  and  viscous  flow  fields.  (2)  The  specific- 
heat  ratios  behind  and  ahead  of  the  wave  front  are 
assumed  constant  and  equal.  This  assumption  is 
not  applicable  to  a  chemically-reacting  flow  where 
the  specific-heat  ratios  for  the  unbumed  and  burned 
gases  are  different.  For  example,  the  specific-heat 
ratio  behind  a  Chapman-Jouguet  (C-J)  detonation  wave 
is  smaller  than  that  of  the  unbumed  gas.  Therefore, 
the  previous  analyses  are  not  applicable  to  actual 


boundary  layers  behind  detonation  waves.  (3)  The 
dynamic  viscosity  coefficient  of  the  gas  behind 
the  wave  front  was  assumed  proportional  to  tem¬ 
perature.  This  assumption  affects  the  wall 
skin-friction  and  heat-transfer  coefficient  but 
not  the  boundary-layer  structure.  (4)  The  Prandtl 
number  is  assumed  constant  across  the  boundary 
layer.  This  assumption  is  reasonable  since  the 
temperature  dependence  of  the  gas  thermal-conduct¬ 
ivity  is  similar  to  the  dynamic-viscosity  coeffi¬ 
cient.  Some  of  these  assumptions  must  be  removed 
and  the  characteristics  of  the  flow  field  must 
be  taken  into  account  for  a  more  realistic  analysis 
of  nonstationary  boundary-layer  flows  behind  deton¬ 
ation  waves.  It  will  be  shown  subsequently  that 
the  available  experimental  heat-transfer  data 
agree  very  well  with  the  present  analysis  for 
planar  flow.  Consequently  this  agreement  also 
adds  confidence  to  the  analyses  of  cylindrical  and 
spherical  detonation  boundary-layer  flows. 

In  this  work,  the  last  three  assumptions  were 
removed.  The  transformed  coordinates  of  Mirels 
and  Hamman  [Ref.  18]  were  modified  in  such  a  way 
that  their  equations  become  applicable  to  flows 
behind  blast  waves  as  well  as  detonation  waves. 

It  has  been  shown  that  direct  initiation  of  a 
detonation  wave  requires  the  instantaneous  release 
of  a  finite  amount  of  ignition  energy  to  guarantee 
a  sufficiently  strong  shock  wave  for  auto-ignition 
in  the  shocked  medium.  When  it  is  ignited  by  an 
exploding  wire,  the  detonation  wave  is,  at  least 
initially,  overdriven,  then  it  decays  rapidly  to  a 
Chapman-Jouguet  (C-J)  detonation  wave  with  a  con¬ 
stant  velocity.  When  such  a  hemispherical  C-J 
detonation  wave  propagates  outward  along  the  major 
diameter,  as  shown  in  Fig.  1,  a  boundary  layer 
forms  behind  it. 

For  C-J  detonation  waves,  similarity  solutions 
exist.  A  complete  description  of  the  method  for 
solving  C-J  detonation  waves  is  given  in  Ref.  13. 

The  detailed  calculations  are  given  by  Saito  [Ref. 

8]  for  detonation  waves  in  2H2+02  mixtures  for 
various  initial  conditions.  For  the  present 
application,  it  was  necessary  to  recalculate  the 
inviscid  flow  functions  and  their  derivatives. 

In  doing  so,  the  compatibility  of  the  solution 
with  the  boundary  conditions  at  the  outer  edge  of 
the  boundary  layer  was  ensured. 

As  pointed  out  in  Refs.  8  and  13  the  detonation 
wave  is  followed  by  a  rarefaction  wave  and  therefore 
all  the  properties  of  the  inviscid  flow  behind  the 
wave  front  decrease  until  the  velocity  approaches 
a  very  small  value  very  near  the  half-distance  that 
the  detonation  wave  has  travelled.  Beyond  that 
point  the  pressure  and  density  remain  constant  and 
a  uniform  stationary  state  results.  However,  in  the 
viscous  flow  behind  a  detonation  wave  the  velocity 
boundary  layer  and  the  thermal  boundary  layer  are 
both  stable  and  grow  with  time  in  the  entire  region, 
except  at  the  origin,  where  the  boundary-layer 
equations  become  singular.  As,  according  to  the 
inviscid  flow,  three  regions  exist  for  the  flow 
behind  a  C-J  detonation  wave,  the  solution  for  the 
boundary- layer  equations  must  be  divided  into  three 
parts  corresponding  to  the  foregoing  three  regions. 
Investigations  were  made,  therefore,  of  boundary 
layers  in  air  behind  strong  blast  waves  and  in 
burned  gases  of  stoichiometric  mixtures  of  hydrogen- 
oxygen  behind  C-J  detonation  waves  for  different 


geometries. 


The  effects  of  Prandtl  number,  the  viscous  ex¬ 
ponent  and  the  wall  surface  temperature  on  the  flow 
structure  are  considered  and  discussed  in  detail. 
Comparisons  have  also  been  made  between  planar  and 
spherical  boundary- layer  growths  behind  blast  waves. 
It  is  shown  that  owing  to  the  rarefaction  wave 
profile  behind  the  shock  fronts  the  densities  drop 
very  rapidly  and  the  boundary  layers  become  very 
thick.  The  spherical  boundary  layer  is  thicker 
owing  to  the  more  rapid  expansion.  This  phenomenon 
does  not  exist  behind  detonation  waves  (or  constant 
speed  shock  waves)  as  the  densities  have  finite 
values. 

This  is  only  a  first  step  in  trying  to  under¬ 
stand  the  viscous  effects  which  exist  in  the  UTIAS 
implosion  facility  at  the  focus  of  an  implosion. 

As  noted  previously,  the  effects  of  viscosity, 
heat  conduction  and  radiation  are  important  for 
obtaining  fusion  in  deuterium.  However,  since 
neutrons  and  y  rays  have  already  been  obtained  in 
this  facility  the  foregoing  effects  cannot  be  too 
limiting. 

Although  the  motivation  was  primarily  for  this 
purpose,  the  present  analytical-numerical  results 
are  general  and  can  be  applied  to  many  boundary- 
layer  flow  problems  behind  blast  and  detonation 
waves. 


2.  1NVISC1P  FLOW  FIELD 
2.1  Basic  Equations 

In  a  mixture  of  chemical  species  of  a  reacting 
gas,  a  great  simplification  can  be  made  in  the 
analysis  of  the  inviscid  flow  by  assuming  momentum 
and  thermal  equilibria  in  the  system.  The  governing 
equations  are  then  reduced  to  three-flow  equations 
for  the  mixture.  The  only  difference  in  the  govern¬ 
ing  equations  between  the  unburned  and  burned  gases 
is  the  equation  of  state.  The  basic  equations  des¬ 
cribing  the  one-dimensional  flow  for  the  mixture 
are. 


If  the  chemical  reaction  is  assumed  to  occur 
within  a  thin  layer  near  the  wave  front,  the  un- 
bumed  and  burned  gases  can  be  treated  as  a  perfect 
gas  with  constant  specific-heat  ratios  r  and  con¬ 
stant  isentropic  exponent  >,  then 

State:  p  =  ;>h  (4) 

The  isentropic  exponent  y  of  the  burnt  2H2*02  is 
about  1.3  at  a  temperature  of  2,000  k  and  1.04  at  a 
temperature  of  20,000  K.  The  pressure  effect  on 
is  small  and  can  be  neglected  when  the  pressure  is 
below  104  atmosphere  [Ref.  24). 

From  similarity  considerations,  it  was  shown 
[Refs. 11, 13]  that  for  a  self-similar  motion  the 
expanding  wave  front  must  either  be  a  power  law 
of  time  or  an  exponential  function  of  time.  For 
most  explosion  problems,  the  power-law  form  is  of 
interest.  Assume  that  the  density  ahead  of  the 
wave  front,  p^,  is  uniform  and  the  wave  front 
position  xs  satisfies 


where  C  and  m  are  constants  and  C  is  determined 
from  the  energy  integral  equation.  The  following 
values  of  m  are  applied  for  explosive  waves: 


m  -  — — — ,  for  a  strong  blast  wave 
o«3  I 

;  (6) 

m  -  1,  for  a  Chapman -Jouguet  wave  | 


Define  the  following  dimensionles  ■  coordinate. 


and  the  dimensionless  variables. 


Mass: 
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Momentum : 
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Energy : 
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where  us  is  the  wave  front  velocity  which  is  a 
function  of  time.  The  following  self-similarity 
equations  result: 


where  a  is  a  numerical  constant  with  values  of  0, 

1,  2  for  planar,  cylindrical  and  spherical  symmetric 
inviscid  flows,  respectively,  p  is  the  mixture  den¬ 
sity,  u  is  the  flow  velocity,  p  is  the  mixture  pres¬ 
sure,  h  is  the  mixture  specific  enthalpy,  t  is  the 
time  and  x  is  the  flow  direction.  The  wall  shear 
stress,  radiation-energy  transfer  and  the  wall  heat 
transfer  are  neglected  for  inviscid  flows. 


[<p-l+S]Re  ♦  ^  '  0  (9) 

[<p-i*q<p^  ♦  =  o  (io) 

♦  yF^  -  2aF  -  yo  ySL  *  o  (11) 
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where  f.  =  0  corresponds  to  the  wave  front  and  (,  =  1 
corresponds  to  the  origin  and  >  is  the  isentropic 
exponent  of  the  gas  behind  the  wave  front.  The  sub¬ 
script  f,  denotes  the  derivative  with  respect  to 
The  parameter  u  is  defined  by 


a 


m-1 

m 


(12) 


where  ct  ~  0  for  a  uniform  wave  front  and  a  <  0  for 
a  decaying  wave  front. 


Initially,  we  assume  that  a  finite  amount  of 
energy  E0  is  released  at  time  zero  in  a  finite 
volume  of  dimension  xQ  in  a  medium  at  rest.  From 
the  energy  integral,  the  constant  C  in  Eq.  (5)  is 
given  by  (Ref .11] 
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where  Kj  is  1,  2i  and  4tt  for  planar,  cylindrical 
and  spherical  symmetry,  respectively,  and  I  is 
defined  by 
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2.2  Boundary  Conditions 
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To  solve  Eqs.  (9)-(ll),  the  boundary  conditions 
immediately  behind  a  shock  or  detonation  wave  front 
have  to  be  determined.  Consider  the  case  where  the 
shock  wave  moves  to  the  right  with  the  velocity  us. 
The  gas  on  the  left  side  of  the  shock,  which  was 
stationary  before  being  traversed  by  the  shock,  is 
given  a  velocity  Uq  to  the  right.  The  situation 
illustrated  above  ma  be  transformed  to  the  case 
where  the  shock  wave  is  stationary.  In  the  new 
frame  of  reference,  he  shock  wave  appears  to  be 
stationary,  while  the  gas  on  the  left  side  of  the 
shock  appears  to  be  flowing  to  the  left  with  the 
velocity  (u0  -  us) .  In  this  frame  of  reference, 
the  conservation  equations  across  the  shock  front 
are  given  by: 


(24) 


For  a  strong  blast  wave  (Ms  °°) ,  the  boundary- 
conditions  at  the  wave  front  (5  =  0)  are 
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For  a  C-J  detonation  wave,  we  have  Ms  *  MfJ  >  ich 
is  defined  by 


Mass : 


p  u  =  o_(u  -u  ) 
00  s  o '  s  o 


(15) 


CJ 


(26) 


Momentum: 


and  s  =  0.  The  boundary  conditions  at  the  wave 
front  are, 


2  2 

p  +pu  =  p  ♦  p  (u  -u  )  (16) 
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Energy : 
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h»  +  7  us2  +  Q  ■  ho  *  1  (us‘u0)2  (17) 

where  Q  is  the  chemical  energy  per  unit  mass  re¬ 
leased  at  the  shock  front.  With  given  initial 
states  (i.e.,  p„,  p„,  h„,  y„  and  Q)  the  state 
immediately  behind  the  shock  front  can  be  expressed 
conveniently  in  terms  of  the  shock  Mach  number  as: 
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2 . 3  Similarity  Solutions 

For  a  strong  blast  wave  the  shock  speed  is 
t ime -dependent ,  hence  the  entropy  change  across 
the  shock  front  decreases  with  time.  The  similar¬ 
ity  conditions  require  an  infinite  strength  wave. 
However,  in  practice  the  strong-shock  approxima¬ 
tion  is  adequate  for  shock  strengths  above  three. 
The  solutions  for  a  strong  blast  wave  can  be 
obtained  by  integrating  Eqs .  (9)-(ll)  with  the 
boundary  conditions  at  the  wave  front,  F.q.  (25). 
Table  1  and  Fig.  2  show  the  dimensionless  flow 
profiles  for  a  strong  spherical  wave  with  T  =  = 

1.4.  Their  corresponding  derivatives  are  shown 
in  Table  2  and  Fig.  3.  The  flow  is  stationary 
only  at  the  origin  (f,  =  1),  which  is  a  singular 
point.  These  similarity  solutions  are  valid 
because  the  pressure  ahead  of  the  blast  wave  can 
be  neglected  compared  with  that  behind  an  intense 
explosion  wave.  The  energy  released  by  the  explo¬ 
sion  and  the  undisturbed  gas  density  play  major 
roles  in  a  strong  blast  wave.  It  is  noted  that 
the  pressure  immediately  behind  the  wave  front 
decreased  inversely  as  its  radius  cubed.  The 
total  energy  contained  in  the  region  behind  the 
blast  wave  is  constant  only  if  p  "t  x‘3.  For 
powers  less  than  or  greater  than  3,  energy  is 
being  added  or  removed  continuously  with  time. 

As  seen  in  Fig.  2,  the  pressure  is  highest  imme¬ 
diately  behind  the  wave  front.  It  then  decays 
as  \  increases,  levels  off,  and  remains  approxi¬ 
mately  constant  for  f,  0.5.  The  particle  velo¬ 
city  varies  in  an  almost  linear  fashion  from  the 
wave  front  to  the  origin  and  vanishes  at  the 
origin.  The  density  approaches  zero  for  £  ■  0.5. 

For  the  C-J  detonation  wave,  it  is  more  con¬ 
venient  to  introduce  a  new  variable  S  defined  as 


Then  the  three  similarity  equations  (9} —(11)  can 
be  reduced  to  two  equations  in  terms  of  <p  and  8 
as  follows  [Ref.  13]: 


For  the  planar  motion  a  equals  zero.  We  can  see 
that  two  solutions  exist:  (1)  (<p-l+£)2-B2  i  0, 
hence  both  and  8 ^  are  zero,  which  corresponds 
to  a  planar  piston  with  constant  speed  driving 
a  constant  velocity  shock  or  detonation  wave. 

(2)  (<f>-lt£)--8“  =  0,  then  both  and  8j  are  finite. 
Integrating  Eqs.  (29)  and  (30)  we  obtain  the 
solution  for  <P(E)  and  6(E)  with  the  proper  signs 
as  follows: 

cp(E)  5  '  ~7T  ’  *  %  (S2) 

e(E)  *  -  *4  £  ♦  e  (33) 


where  (for  the  C-J  condition) 

90  *  S0  =  1  (34 

The  subscript  o  denotes  the  value  at  f  =  0.  Using 
Eqs.  (28)  and  (31),  we  get 


F  =  ~  R  (36) 

If  ?0  and  d0  are  known,  using  Eqs.  (32),  (33),  (35) 
and  (36),  cp,  R,  F  can  be  calculated;  <j^  and  S0  are 
determined  from  Eqs.  (27)  and  (28). 

For  the  cylindrical  and  spherical  detonation 
wave,  5  4  0,  *  80  =  1 .  This  results  in  9r  -*■  -°° 

and  fir  ■»  -<»  at  the  front,  £  =  0.  Therefore,  we 
cannot  integrate  Eqs.  (29)  and  (30)  directly  from 
the  boundary  of  the  C-J  detonation  wave  front.  To 
overcome  this  difficulty,  the  perturbation  method 
is  used  [Ref. 13].  The  solution,  in  the  form  of  a 
perturbation  series,  in  the  immediate  neighbourhood 
of  the  C-J  front  is  given  by  the  following  expres¬ 
sions  (with  proper  signs): 
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In  the  C-J  detonation  wave,  the  shock  front  propa¬ 
gates  at  a  constant  velocity  and  the  entropy 
increase  is  identical  for  all  gas  particles  across 
the  wave  front .  The  invariant  then  reduces  to  the 
following  form: 


-=-  «  K,  »  constant 
F  2 


Kg  can  be  evaluated  using  the  boundary  conditions  at 
the  front,  £  «  0. 


With  Eqs.  (37)  and  (38)  the  values  <p(o+A£)  and 
8(o+A£)  can  be  calculated  from  the  values  <p0  and  60, 
where  AC  is  very  small.  From  these  values,  q>(o*A£) 
and  8(o+A£),  Eqs.  (29)  and  (30)  can  then  be  inte¬ 
grated  to  obtain  <J>(£)  and  6(E).  Using  Eqs.  (35) 
and  (36),  R(E)  and  F(E)  are  also  obtained. 

A  stoichiometric  mixture  of  hydrogen  and  oxygen 
was  used.  The  specific  heats,  specific-heat  ratios, 
isentropic  exponent  and  equilibrium  sound  speed  for 
the  unbumed  and  burned  gases  are  given  by  Benoit 
[Ref.  24], 

In  the  Chapman- Jouguet  detonation  wave  the  shock 
front  propagates  at  a  constant  velocity  and  the 
entropy  increase  for  all  gas  particles  across  the 


wave  front  are  identical.  Table  3  and  figs.  4-6 
show  the  dimensionless  pressure,  velocity  and 
density  profiles  behind  the  planar,  cylindrical 
and  spherical  C-J  detonation  waves  with  initial 
conditions  of  =  13.6  atm,  =  1.4  and  I  =  300  K 
The  corresponding  derivatives  which  are  needed  in 
the  boundary-layer  analysis  are  shown  in  Table  4 
and  Figs.  7-9.  It  is  noted  that  there  are  three 
regions.  From  the  wave  front  (E,  =  0)  to  f.  =  0.48 
is  the  expansion  region  and  from  f.  =  0.48  to  the 
origin  is  the  stationary  tone.  The  particle 
velocity  decreases  sharply  from  its  maximum  value 
at  f,  =  0  to  zero  at  i.  =  0.48.  Similar  to  the  blast 
wave,  the  pressure  behind  the  wave  front  decreases 
as  f,  increases  and  then  remains  constant  for 
E  '  0.48.  However,  the  gas  density  resembles  the 
pressure  profile  and  does  not  approach  zero  close 
to  the  wave  midpoint.  It  can  be  seen  from  Figs. 

7-9  that  a  discontinuity  exists  in  the  flow  deriva¬ 
tives  at  the  position  where  the  particle  velocity 
is  zero.  The  flow  derivatives  also  approach  an 
infinite  value  at  the  detonation-wave  front.  As 
mentioned  previously,  at  the  origin  the  basic 
equations  are  singular. 


the  x  and  y-dirc  tions,  respectively,  and  Fr  is  the 
Prandtl  number.  The  temperature  dependence  of  the 
dynamic-viscosity  coefficient  u  is  assumed  to  be  of 
the  form: 


where  the  subscript  r  denotes  a  reference  value  and 
a  is  the  viscous  exponent. 

3 . 2  Transformed  Boundary- Layer  hquat ions 

The  following  transformed  coordinates  were  used 
by  Mirels  et  al . 
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3.  VISCOUS  FLOW  F1EI.D 

3.1  Governing  Hquat ions  for  laminar  Boundary -layer 
Flows  — 


where  A  =  2mF„C‘’'°*^v  -_/p  ,  v  =  u  /S  .  As  men 
tioned  previously  this  set  of  transformed  coordi- 
nates  cannot  be  used  in  chemical -reacting  flows. 
By  defining  the  transformed  coordinate  n  as 


The  laminar  boundary  layer  behind  a  strong 
blast  wave  moving  with  nonuniform  velocity  was 
investigated  by  Mirels  et  al  [Refs. 18,  23]  and 
Chen  and  Chang  [Ref.  19] .  The  laminar  boundary 
layer  behind  a  detonation  wave  has  not  been  solved 
as  yet.  As  mentioned  in  the  Introduction,  the 
last  three  assumptions  made  by  Mirels  et  al  [Ref. 

18]  are  not  applicable  to  a  detonation -wave  induced 
boundary  layer.  In  order  to  generalize  the  analysis, 
these  assumptions  are  modified:  (1)  The  specific 
heats  and  the  specific-heat  ratios  behind  and  ahead 
of  the  wave  front  are  different ,  but  constant  in 
each  region.  (2)  The  Prandtl  number  may  be  variable 
across  the  boundary  layer.  (3)  The  temperature 
dependence  of  the  dynamic-viscosity  coefficient 
of  the  gas  behind  the  wave  front  has  a  power  law 
variation  with  temperature.  Other  assumptions  are 
made  that  the  chemical  reaction  takes  place  only  at 
the  C-J  detonation-wave  front  and  that  the  gases 
behind  and  ahead  of  the  detonation  wave  are  perfect. 
Therefore,  the  equations  governing  the  laminar 
boundary-layer  flow  induced  by  a  wave  front  are  as 
follows: 
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where  A  =  2C2  'm‘J  'p^bF^/p  C  is  given 

by  Eq.  (13),  F0  is  given  by  Eqs.  (25)  or  (27)  and 
b  =  l7h„-l)/T„(Y  -l)]Cp„/Cp  is  the  molecular  weight 
ratio  of  the  gases  behind  and  ahead  of  the  wave 
front,  then  n  is  applicable  to  chemical -reacting 
flows  as  well  as  non-reacting  flows.  It  can  be 
shown  that  when  u>  *  1 ,  b  =  1  and  y  *  then  A  and 
r,  become  identical  to  those  given  by  Mirels  and 
Hamman  [Ref. 18]. 


Define  the  following  dimensionless  variables: 


m,  n) 


g(E.  n)  =  p- 
e 


(46) 


where  the  subscript  e  denotes  the  edge  of  the  boun¬ 
dary  layer.  Since  the  pressure  is  constant  in  the 
y-direction,  the  density  ratio  is  given  by 


and  o  is  a  boundary-layer-type  parameter,  with  a 
value  0  for  plane  flow  and  1  for  axisymmetric 
flow.  The  pressure  p  is  constant  in  the  y-direc¬ 
tion,  which  is  normal  to  the  surface  or  flow- 
direction  x.  The  velocities  u  and  v  are  along 
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If  a  scalar  stream  function  ip  is  now  introduced 
such  that 
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fhc  continuity  equation  is  then  satisfied  automatic 
cally.  It  is  not  difficult  to  show  that  the  rela¬ 
tionship  between  the  scalar  stream  function  i|>  and 
the  dimensionless  variable  f(f.,  n)  is  given  by 
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Applying  the  transformed  coordinates,  dimension¬ 
less  variables  and  the  relationship  between  and 
f  to  l:.qs.  (39)  to  (41),  the  transformed  boundary- 
layer  equations  are  then  obtained, 


where  the  subscript  w  denotes  the  wall  surface.  It 
should  be  noted  that  from  the  definition  of  gw 
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should  be  a  function  of  f  only,  for  a  self-similar 
boundary -layer  flow.  Therefore,  two  cases  exist: 
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(1)  For  a  strong  blast  wave,  according  to  its 
definition,  we  have 
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where  a  is  defined  by  Eq.  (12),  the  values  F,  R, 
ip  and  their  derivatives  Fjr ,  and  are  obtained 
from  Eqs.  (9)  to  (11),  the  subscripts  n  and  £ 
denote  the  derivatives  with  respect  to  n  and  5, 
and  B  is  defined  by 
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It  can  be  shown  that  if  u)  =  1  (i.e.,  B  =  1)  and 
Pr  *  constant,  Eqs.  (SO)  and  (51)  reduce  to  those 
derived  by  Mirels  and  Hamman  [Ref. 18).  (The 
details  can  be  found  in  Appendix  A.) 

The  required  boundary  conditions  are  given  by 
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After  the  transformation,  they  become 


W.  0)  -  fnU,  0)  =  0 

g(£.  0)  =  gw(C)  (55) 


so  that  gw  =  0 . 

2 

(2)  For  the  C-J  detonation  wave,  we  have  Ms^  =  Mqj  = 
constant.  Consequently,  gw  is  a  function  of 
£  only. 

It  should  be  noted  that  for  the  strong  blast  wave 
or  for  the  so-called  cold-wall  model  (i.e.,  Tw  =  0) 
the  viscous  exponent  u  has  to  be  unity.  Otherwise, 
Eqs.  (50)  and  (51)  will  be  singular  at  the  wall. 

3.5  Numerical  Method  and  Procedure 


The  nonlinear  equations,  (50)  and  (51),  with 
the  boundary  conditions  given  by  Eqs.  (55)  and  (56), 
were  solved  numerically  by  an  implicit  six-point 
finite-difference  scheme  developed  in  Ref.  25.  The 
initial  profiles  are  required  for  a  finite-differ¬ 
ence  method.  Consequently,  the  solution  at  the 
start  of  the  boundary  layer  has  to  be  obtained 
first.  At  C  »  0,  the  partial-differential  equa¬ 
tions,  Eqs.  (50)  and  (51),  become  a  set  of  ordinary- 
differential  equations  as  follows: 


4.  •*?£<] 

♦  (n-«Pf)gr|  *  o 


(58) 


(59) 


Equations  (58)  and  (59)  were  solved  using  the 
Newton-Raphson  method,  or  a  shooting  method. 

Equations  (50)  and  (51)  are  linearited  in  a 
form  suitable  for  an  iteration  scheme  by  introducing 
the  following  function 


W 


3f 

3n 


(60) 


It  is  really  only  required  for  the  momentum  equa- 


n  =  0: 


tion  and  not  for  the  energy  equation.  "fie  momentum 
equation  in  terms  of  W  can  be  written  as 


(1-C)2°R1'u,(F/Fo)Bw‘P) 


*  |(l-?;)20R1'“(F/F0)u‘Bn  ♦  (n-qjf) 
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and  the  energy  equation  in  the  linearized  form  is 


U-^V-CF/F^  ~  g‘P> 


layer  are  not  satisfied.  In  order  to  avoid  a  dis¬ 
continuity  in  the  gradients  of  the  dependent 
variables  at  the  edge  of  the  boundary  layer,  they 
were  forced  to  select  a  curve  for  ne  in  their  num¬ 
erical  procedure.  The  boundary- layer  profiles  were 
then  recalculated  using  the  selected  ne  distribution 
However,  in  the  present  work  the  inviscid  flow 
equations,  Eqs.  (9) -(11),  and  the  boundary- layer 
equations,  Eqs.  (50)  and  (51),  are  solved  simul¬ 
taneously  with  the  same  step  size  A£.  No  inter¬ 
polation  errors  are  introduced  and  the  compatibility 
conditions  are  satisfied.  The  selection  of  an  ne 
distribution  is  not  necessary  in  this  work.  The 
advantage  of  the  present  numerical  procedure  is  the 
reduction  of  the  number  of  convergent  iterations 
for  the  nonlinear  differential  equations  and  the 
consequent  computation  costs. 

In  the  calculations,  BO  mesh  points  were  used 
in  the  n-direction  with  X  =  1.02  and  an  initial  step 
size  An  =  0.026.  The  step  size  AC  =  0.002  was 
applied.  Normally  two  or  three  iterations  were 
required  to  converge  the  solution  within  an  error 
of  10-6. 


*  (l-C)2aR1'u(F/Fo)ul  (  j?-  )n  ♦  (n-<pf) 

-  2C  [f^ff-  -  ^  [2o*a(2w-l)l]  g<p) 

-2«  [f*  (1'^)  (S-Tr)]s(P) 

-  2C(l-C-<pw)g£p3  -  (1-C)2°  ~  R1‘<“(F/F0)‘W 

(62) 

where  superscript  p  denotes  the  order  of  the  itera 
tion  process  and  the  quantities  without  the  super¬ 
script  denote  those  evaluated  at  the  (P-1)  itera¬ 
tion  order. 

A  computer  program  was  written  based  on  the 
numerical  scheme  outlined  in  Ref.  25  and  Eqs.  (61) 
and  (62)  were  solved  for  four  cases  of  practical 
interest.  The  results  are  presented  in  Section  4. 
In  this  method,  either  equal  intervals  or  nonequal 
intervals  can  be  used  in  the  n-direction.  The 
interval  in  the  n-direction  is  increased  in  a  geo¬ 
metric  progression  as 


where  X  is  a  constant  set  with  a  value  slightly 
greater  than  unity  and  i  is  the  index  of  the  n- 
coordinate  for  the  difference  net.  There  are  two 
differences  in  the  numerical  procedure  between 
the  present  work  and  Liu  and  Mirels  [Ref.  23] . 
First,  in  this  work  the  momentum  equation  (SO)  is 
transformed  into  a  second-order  differential  equa¬ 
tion  through  a  function  defined  by  Eq.  (60).  The 
numerical  method  of  Blottner  [Ref.  26] ,  which  was 
applied  by  them,  introduced  a  transformed  normal 
velocity  and  retained  the  continuity  equation  in 
order  to  avoid  third-order  derivatives  in  the 
momentum  equation.  Second,  Liu  and  Mirels  used 
an  interpolation  method  in  the  inviscid  flow 
values  which  appear  in  the  coefficients  of  the 
boundary- layer  equations.  Consequently,  the  com¬ 
patibility  conditions  at  the  edge  of  the  boundary 


3.4  Boundary-Layer  Characteristics 

The  relations  between  the  transformed  coor¬ 
dinates  (£,  n)  and  the  physical  coordinates  (x,y,t) 
are  given  by  Eqs.  (43)  and  (45), 


where 
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The  boundary-layer  thickness  6e,  displacement 
thickness  6*  and  the  momentum  thickness  6  are 
expressed  by 


0 

f°°  PU 

J 


where 


*  H,H2  |  6  gdn 

(67) 

0 

Peue  J  121 

(68) 

I1  -  IT]  *H1H2S2 

l  e  > 

*  {  (R-fn)dn  (69) 

o 

A 

■  j  yi-ydn  (70) 

o 


7 


and  rie  is  the  value  of  n  at  the  edge  of  the  boun¬ 
dary  layer. 


Ref.  30: 


From  the  definition  of  the  shear  stress  and 
the  skin-friction  coefficient  Cf: 
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the  dimensionless  skin  friction  coefficient  is 
expressed  by 
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where  the  Reynolds  number  Re  is  given  by 
P  *  ^ 


Re  = 
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The  heat  transfer  at  the  wall  is  given  by 
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i  =  1 


for  the  thermal  conductivity  of  the  burned  gas,  and 
n 

-  N.p. 

(79) 
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for  the  viscosity  of  burned  gas,  where 
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and  N,  u,  K,  M  are  the  mole  fraction,  viscosity, 
conductivity  and  molecular  weight  of  the  species, 
respectively.  Subscripts  i  and  j  denote  the  ith 
and  jth  species.  The  details  of  the  calculation 
are  given  in  Appendix  B. 

3 .6  Discussion  of  Boundary-Layer  Equations 


where  K  is  the  thermal  conductivity  of  gas. 
the  definition  of  Stanton  number  St: 


With 
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where  Hg  =  he  ♦  (l/2)ue2-  The  normalized  expression 
for  wall  heat  transfer  is 


As  stated  before,  three  regions  are  considered 
in  treating  the  boundary  layer  behind  a  C-J  detona¬ 
tion  wave.  For  the  expansion  region  from  £  »  0  to 
£  v  0.5,  the  inviscid  flow  plays  a  significant  role 
in  the  development  of  the  boundary  layer.  Equa¬ 
tions  (50)  and  (51)  are  applied  directly.  For 
the  stagnation  region  from  £  v  0.5  to  £  *  0.95, 
because  the  velocity  approaches  a  very  small  value, 
Eqs.  (50)  and  (51)  become 
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3.5  Properties  of  Hydrogen-Oxygen  Stoichiometric 
Mixtures 


A  stoichiometric  mixture  of  hydrogen  and  oxygen 
is  taken  as  the  working  fluid.  The  properties  of 
the  unburned  and  burned  gases,  i.e.,  the  gases 
behind  and  ahead  of  the  detonation  wave,  are  taken 
from  Ref.  24.  The  chemical-reaction  equation  is 
as  follows: 


(i-O20r1'u‘  (  jr  )  (  &  8n)n  *  % 

2C  {-  £  [2o*a(2u>-l)]gn  ♦  &  g  ♦  (1-OgJ  ...  (82) 


In  the  region  from  £  *  0.95  to  £  *  1,  the  solution 
is  given  by 


2H2*02  7  NjH20  +  N20H  ♦  N3H2  *  Nj02  ♦  NjH  ♦  Nfi0 

where  Nj,  Nj,  ....  denote  molar  concentrations  of 
HgO,  OH,  ....  in  the  burned  gases. 

The  initial  conditions  of  the  stoichiometric 
mixture  of  hydrogen  and  oxygen  are  p„  -  13,6  atm, 

T„  »  298. 15°K,  y»  *  1.4  and  us  >  2.9815x10s  cm/sec. 

The  thermal  conductivity  and  kinetic  viscosity 
for  the  unbumed  and  burned  gases  are  calculated 
using  the  semi-empirical  relations  outlined  in 


T(y,t)  -  (Te-Tw)erf(y/^Dt)  ♦  T,  (83) 

where  T  is  temperature,  erf  (y/*'4Dt)  is  the  error 
function  and 
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where  C_  is  the  specific  heat  of  the  gas.  The  tem¬ 
perature  gradient  at  the  wall  decreases  in  propor¬ 
tion  to  f'/2.  in  the  initial  phase,  the  tempera¬ 
ture  gradient  is  so  large  that  a  temperature  jump 
exists  at  the  wall.  Details  of  the  derivation  of 
Eq.  (83)  are  given  in  Appendix  C. 

Three  interesting  observations  may  be  made 
regarding  Eqs.  (50)  and  (51).  First,  as  mentioned 
above,  the  equations  are  singular  at  (  =  1.  The 
solution  near  £  =  1  is  more  of  analytical  interest, 
as  the  boundary-layer  concept  is  not  valid.  Second, 
the  effect  of  the  viscous  exponent  ui  on  the  boun¬ 
dary  layer  depends  on  the  inviscid  solutions  R 
and  (F/F0) .  When  u)  <  1 ,  the  gas  viscosity  and 
thermal -conduct ion  effects  increase  as  w  increases. 
Third,  the  coupling  between  the  momentum  and 
energy  equations  is  through  B  and  the  last  term  on 
the  right-hand  side  of  Eq.  (50)  and  the  second 
term  on  the  right-hand  side  of  Eq.  (51). 

It  is  convenient  to  rewrite  Eqs.  (50)  and 
(51)  in  the  general  form 


a,W  „  +  a,K  ♦  a,W  +  a, 
l  nn  2  n  3  4 


a5W£ 
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where  Che  are  the  coefficients  and  W  stand  for 
fr)  and  g.  Equation  (85)  has  the  same  form  as  the 
heat  conduction  equation.  The  coefficient  as, 
which  plays  the  same  role  as  the  thermal  diffus- 
ivity  in  the  heat -conduct ion  equation,  is  negative 
in  regions  of  reverse  flow.  For  axisymmetric 
boundary-layer  flows  (i.e.,  a  *  1),  a[  becomes 
very  small  when  £  =  1 .  It  means  that  the  viscosity 
and  thermal  conduction  effects  become  very  small 
when  tq  s  0.  When  aj  =0,  Eq.  (85)  is  no  longer 
parabolic,  but  hyperbolic.  Consequently,  all 
numerical  methods  for  solving  the  boundary- layer 
equations  fail  in  this  case.  The  simple  heat- 
conduction  equation  can  then  be  applied. 

The  unsteady  effect  on  the  boundary- layer 
development  can  be  estimated  as  follows.  Consider 
the  boundary  layer  induced  by  a  uniform  wave  (i.e.  , 
wave  velocity  is  constant,  or  «  =  1).  The  steady- 
state  boundary- layer  equations  are: 


(1-£)2°R1~W  I  ]  (Bf  ) 


nn  n 

«P  f 

£  -  n£ 


!§* 

Rq> 


(1-S)2cR 


-(£)[( 


Pr  gnln 


J5SL  Bf2 1 
y  F  nnj 


h  fS  \ 
g  yF  r 


(86) 


'  f«n 


g  (87) 


Comparing  Eqs.  (50)  and  (51)  with  m  *  1  and  Eqs. 
(86)  and  (87),  it  can  be  shown  that  the  unsteady 
effect  can  be  neglected  only  when  the  following 
conditions  are  satisfied: 


(1)  (l*2£o)n  « 


(2)  l-£  «  <pf 

T  n 

The  nonuniformity  of  the  wave  front  increases 
the  unsteady  effect  on  the  boundary -layer  structure. 
The  order  of  magnitude  for  the  unsteady  terms  depends 
on  the  inviscid  flow  characteristics  and  it  is  dif¬ 
ficult  to  estimate. 

4.  BOUNDARY -LAYER  DEVELOPMENT  BEHIND  EXPLOSIVE 
WAVES 

4 . 1  Boundary  Layer  Behind  Strong  Blast  Waves 

Three  sets  of  solutions  were  obtained.  The 
first  set  was  for  the  boundary  layer  induced  by  a 
spherical  strong  blast  wave.  In  order  to  check 
the  numerical  procedure  applied  here,  the  case  with 
y  =  'Y0o=  1.4,  ui  ~  1,  b  =  1,  g w  =  0  and  Pr  =  0.72 
behind  a  spherical  blast  wave  was  examined.  Re¬ 
sults  were  obtained  beyond  £  =  0.9S.  Figures  10 
and  11  show  the  normalized  velocity  and  temperature 
profiles.  For  comparison,  the  corresponding  results 
obtained  by  Liu  and  Mirels  (Ref.  23]  are  shown  in 
Fig.  12.  It  can  be  seen  that  agreement  between 
the  two  sets  of  results  is  excellent.  Further 
tabular  comparisons  are  given  in  Tables  5  and  6. 
Agreement  of  the  present  results  with  the  series 
expansion  results  (Ref. 18]  for  the  wall  derivatives 
at  £  =  0  is  better  than  those  obtained  by  Liu  and 
Mirels  [Ref.  23]. 

The  thermal  boundary- layer  thickness  is  larger 
than  the  velocity  boundary- layer  thickness  since 
Pr  s  i  was  applied.  As  £  increases,  the  velocity 
gradient  increases  at  a  rate  which  is  larger  than 
that  for  the  temperature  gradient.  Consequently,  a 
velocity  overshoot  occurs.  Results  of  other  cases 
with  different  Pr  indicate  that  (1)  the  velocity 
overshoot  decreases  as  Pr  increases,  (2)  there  is 
no  flow  reversal  when  the  Pr  is  varied  from  0.72 
to  more  than  3,  because  the  density  decreases 
rapidly  with  £.  It  can  be  seen  from  Figs.  10  and 
11  that  the  velocity  and  thermal  boundary- layer 
thicknesses  decrease  with  He  as  £  increases  (the 
actual  thickness  increases  with  £) .  They  become 
zero  at  the  origin  (£  =  1).  Therefore,  fn  =  1 
for  all  n  at  £  -  1 . 

4 . 2  Boundary  Layer  Behind  C-J  Detonation  Waves 

A  second  set  of  solutions  was  determined  for 
the  C-J  detonation  wave  in  a  stoichiometric  mixture 
of  hydrogen-oxygen  with  ic  =  0.75,  Pr  =  2.26,  Tw  * 

T  =’300%  Y  -  1.14,  Y  *  1.4  and  p  =  13.6  atm 
for  different  flow  geometries.  The  distance  used 
from  the  origin  to  the  wave  front  is  10  cm  (corre¬ 
sponding  to  the  implosion  chamber  wall)  and  the 
wave  velocity  is  2982  m/s.  Figures  13  to  20  show 
the  normalized  velocity  and  temperature  distribu¬ 
tions  across  the  boundary  layer  for  the  planar, 
cylindrical  and  spherical  inviscid  flows.  The  wall 
derivatives  are  shown  in  Table  7  and  in  Figs.  21  to 
24.  From  the  results  presented  above,  the  following 
observations  are  made:  (1)  Flow  reversal  occurs  for 
all  cases  with  Fr  *  2.26  when  £  is  greater  than  0.4. 
The  magnitude  of  the  flow  reversal  increases  as  o 
and  o  increase.  As  a  result,  it  has  a  maximum  flow 
reversal  for  the  spherical  flow  geometry,  i.e.,  for 
5=2  and  o  *  1.  (2)  The  velocity  gradient  at  the 

wall  increases  as  £  increases  and  reaches  a  maximum 
value  at  £  between  0.2  and  0.25  and  then  decreases 
with  £  increasing.  (3)  For  the  case  with  o  «  1, 
the  wall  temperature  gradient  increases  at  a  rate 


greater  than  for  o  =  0,  as  E  increases.  (4)  The 
thickness  of  the  velocity  boundary  layer  is  greater 
than  that  of  the  thermal  boundary  layer. 

A  third  set  of  solutions  was  obtained  for  a 
spherical  C-J  detonation  wave  in  a  stoichiometric 
mixture  of  hydrogen -oxygen  (T^  =  300°K,  >  =  1.14, 
y  =1.4  and  p  =  13.6  atm)  with  different  thermal 
propert ies . 

As  mentioned  before,  for  Pr  %  1 ,  flow  reversal 
exists  when  E  '  0.3.  In  order  to  examine  the 
effect  of  the  Prandtl  number  on  the  flow  structure, 
other  cases  with  Pr  =  1.5,  Pr  =  1.0  and  Pr  =  0.72 
were  also  studied.  Figures  25  to  30  show  the 
normalized  velocity  and  temperature  for  different 
Prandtl  numbers.  From  these  figures,  including 
Fig.  19,  it  can  be  seen  that  for  Pr  '  1 ,  the  mag¬ 
nitude  of  the  flow  reversal  decreases  with  Pr 
decreasing,  whereas  for  Pr  <  1 ,  the  flow  reversal 
disappears  and  a  velocity  overshoot  occurs.  The 
magnitude  of  velocity  overshoot  increases  as  the 
Pr  decreases.  It  is  clearly  shown  in  Fig.  27  that 
for  Pr  =  1  neither  a  flow  reversal  nor  a  velocity 
overshoot  exists.  Figures  31  and  32  show  the  wall 
derivatives  of  velocity  and  temperature,  respec¬ 
tively.  It  can  be  seen  from  Fig.  31  that  when 
Pr  <  1 ,  fnn(f.,0)  increases  monotonically  with  E 
increasing.  However,  when  Pr  ■  1,  fnq(f.,0) 
increases  at  first  and  reaches  a  maximum  value 
at  about  E  =  0.23,  then  decreases  to  a  negative 
value.  Separation  is  said  to  occur  when  the 
derivative  fnn(E,0)  changes  from  a  positive  to  a 
negative  value.  The  boundary- layer  thicknesses 
of  velocity  and  temperature  are  shown  in  Fig.  33. 

It  is  shown  that  when  Pr  >  1 ,  the  velocity  boundary 
layer  thickness  increases  slightly  as  Pr  increases, 
whereas  the  thermal  boundary- layer  thickness  de¬ 
creases  considerably.  It  implies  that  when  Pr  >  1 , 
more  kinetic  energy  of  flow  is  lost  from  the  main¬ 
stream  than  when  Pr  <  1 .  Due  to  imbalance  between 
the  velocity  and  thermal  boundary- layer  thicknesses, 
the  flow  is  forced  to  separate  at  the  wall  or  to 
have  a  velocity  overshoot  at  the  edge  of  the  boun¬ 
dary  layer.  This  separation  is  not  observed  in 
the  blast-wave  case. 

Table  8  shows  the  variation  of  the  thicknesses 
of  velocity  and  temperature  with  Pr.  It  was  found 
that  the  ratio  of  the  velocity  boundary-layer  thick¬ 
ness  to  the  thermal  boundary- layer  thickness  is 
proportional  to  the  square  root  of  the  Prandtl 
number,  as  shown  for  a  steady  plane  boundary-layer 
flow  [Ref.  27],  For  the  spherical  C-J  detonation 
wave,  it  can  be  seen  from  Table  8  that  when  Pr  >  1 , 
the  velocity  boundary-layer  thickness  is  greater 
than  for  the  temperature  layer.  When  Pr  <  1,  the 
converse  statement  is  true  and  when  Pr  =  1,  the 
two  thicknesses  are  identical. 

The  effect  of  the  viscous  exponent  and  the 
surface  temperature  on  the  boundary-layer  flow  was 
also  investigated.  The  normalized  velocity  and 
temperature  profiles  are  shown  in  Figs.  34  and  35 
for  the  case  of  u)  *  1,  and  in  Figs.  36  and  37  for 
the  case  of  Tw  «  0,  respectively.  The  effects  of 
<e  and  T„  on  the  velocity  and  temperature  gradients 
are  shown  in  Table  9  and  in  Figs.  38  to  41.  It 
can  be  seen  from  Table  9  and  Figs.  38  to  41  that: 

(1)  when  uj  changes  from  0.75  to  1,  the  velocity 
and  temperature  gradients  on  the  wall  increase  by 
a  factor  of  nearly  1.8.  (2)  The  wall-temperature 

effect  on  the  wall  gradients  is  not  significant. 


The  thermal  boundary- layer  thickness  growing 
on  the  wall  of  the  major  diameter  of  the  UT1AS 
hemispherical  implosion  chamber  was  obtained  numer¬ 
ically  from  Eqs.  (50)  and  (51),  and  analytically 
from  Eq.  (83).  Figure  42  shows  the  variation  of 
the  thermal  boundary- layer  thickness  with  E  just 
when  the  detonation  wave  reaches  the  hemispherical 
wall  of  the  implosion  chamber  at  34  usee  from  init¬ 
iation,  and  Fig.  43  shows  the  variation  with  time 
at  the  origin  for  the  case  of  T„,  =  300°K,  b  =  1.26, 

F  =  0.1511,  R  =  0.67,  US  =  2982  m/s  and  xs  =  10  cm. 
The  predicted  thermal  boundary- layer  thickness 
using  the  analytical  solution  is  about  0.054  mm  at 
the  origin,  which  is  very  close  to  0.058  mm  predicted 
numerically. 


4 . 3  Some  Discussions  on  Boundary-Layer  Thickness 


Figure  44  shows  the  variation  of  boundary- layer 
thicknesses  for  C-J  spherical  and  planar  detonation 
waves  in  2112+02,  and  strong  spherical  and  planar 
blast  waves  in  air  under  the  same  initial  conditions 
and  wave  velocity  us  and  wave  position  xs.  For  a 
comparison,  the  constant  speed  nonstationary  planar 
shock-wave  case  in  air  is  also  shown  in  Fig.  44. 

The  initial  conditions  are:  p„  =  13.6  atm,  T,*,  = 
300°K,  wave  speed  us  =  2981.7  m/sec,  and  wave 
position  xs  =  10  cm.  It  can  be  seen  from  Fig.  44 
that  the  variation  of  the  boundary- layer  thickness 
with  distance  behind  the  wave  behaves  quite  differ¬ 
ently  for  a  C-J  detonation  wave,  a  strong  blast  wave 
and  a  constant  speed  planar  shock  wave.  It  is  help¬ 
ful  to  recall  the  relation  for  the  boundary-layer 
thickness,  Eq.  (67), 


<*„  I"  [1/2 

b  —  F 
Poo  J 


(2£.) 


1/2 


(1-C)°R 


ne 

j  gdn  (67) 

0 


in  order  to  understand  these  variations.  We  can  see 
from  Eq.  (67)  that  the  density  ratio  plays  a  very 
important  role  in  the  development  of  the  boundary 
layer  behind  a  given  wave  front.  For  a  C-J  detona¬ 
tion  wave  when  E  >  0.5,  the  velocity  of  the  inviscid 
flow  becomes  very  small,  but  it  does  not  matter 
because  the  density  and  pressure  remain  at  reason¬ 
able  levels.  For  a  constant  speed  planar  shock 
wave,  a  «  0,  R  =  constant  and  ^)r'e  gdn  is  also  con¬ 
stant  (that  is,  similarity  exists  and  the  boundary- 
layer  profiles  are  identical  anywhere),  therefore 
the  boundary-layer  thickness  S  can  be  expressed 
as 


6  =  c/E 

where  c  is  a  constant.  For  the  strong  blast  wave, 
however,  when  E  >  0.5,  the  density  ratio  becomes 
very  small  (as  shown  in  Table  1  and  Fig.  2).  It 
means  that  when  E  *  0.5,  the  inviscid  flow  is  so 
rarefied  that  the  boundary  layer  gets  comparatively 
thick  (as  shown  in  Fig.  44).  Consequently,  the 
assumptions  made  in  derivation  of  the  boundary- 
layer  equations  would  cease  to  apply  and  the  plots 
in  Fig.  44  for  these  cases  for  larger  distances  from 
the  wave  are  no  longer  of  physical  interest. 


4.4  Heat  Transfer  to  the  Wall  behind  a  Detonation 
Wave  in  2H,+0,  Mixture 

The  relation  for  calculating  wall  heat  transfer 


10 


behind  a  C-J  detonation  wave  moving  into  a  station-  Ref-  27  over  the  entire  200  usee .  However,  it 

ary  mixture  of  2H,«0,  is  given  by  (see  Appendix  D)  only  agrees  well  with  the  data  of  Ref.  15  from 

about  ISO  to  200  usee.  For  structure-independent 
_(j#5  detonations,  as  shown  in  Fig.  49,  velocity,  pressure 

%  =  *®i'®2^t)  (68)  an<j  density  behind  the  wave  front  decrease  rapidly, 

then  the  values  of  p  and  p  level  off  to  a  finite 

where  value  at  about  £,  •  0.5.  However,  the  velocity 

after  C,  -  0.5  becomes  so  small  that  only  conduction 
1/2  heat  transfer  becomes  dominant.  That  may  be  the 

g  =  I  w  n  w  II  l  (l-f,)°  (69)  main  reason  why  the  present  laminar  boundary- layer 

*  7"^  /2  pr  '  l  F0*  J  J  analysis  gives  a  good  representation  of  the  varia¬ 

tion  of  heat  transfer  to  the  wall  with  time  behind 
C-J  detonation  waves. 

B2  =  uu(p»KJ°'5  (7°) 

5.  CONCLUSIONS 

It  can  be  seen  that  B(  is  a  function  of  f)  only,  it 

decreases  as  £  increases  from  0  to  1 .  For  planar  A  modification  to  the  transformed  coordinates 

flow  with  o  =  0  it  reduces  to  the  following  form  of  Mirels  and  Haraman  [Ref.  18)  makes  it  possible 

to  apply  the  transformed  equations  to  flows  behind 
gu~l(g  )  .  .  ltoil/2  either  blast  or  Chapman-Jouguet  detonation  waves, 

g  =  1  w  n  w  IX  l  (71)  It  was  shown  that  the  Prandtl  number  controls  a 

1  Y'l  pr  '  v  R0R  J  J  reversed-f low  phenomenon  at  the  wall  and  a  velocity 

overshoot  near  the  edge  of  the  boundary  layer.  No 

while  B2  represents  the  effects  of  the  initial  flow  reversal  occurs  in  the  blast -wave  case.  The 

conditions.  For  the  structure-independent  detona-  effect  of  Prandtl  number  on  flow  reversal  in  the 

tion  U(u  is  a  function  of  the  initial  conditions  blast-wave  case  is  not  as  significant  as  in  the 

only.  From  the  expression  for  Bi,  it  can  be  seen  detonation-wave  case.  The  magnitude  of  the  velocity 

that  qw  is  also  a  function  of  position.  For  a  overshoot  for  the  Chapman-Jouguet  detonation-wave 

given  position  x,  5  can  be  expressed  as  case  is  larger  than  that  for  the  blast -wave  case. 

The  viscous  exponent  u  has  a  significant  effect  on 
u  t  the  wall  skin  friction  and  heat  transfer.  For 

£  *  _ ^2)  w  <  1,  the  value  of  wall  derivatives  reduces  by  a 

x*uCjt  factor  of  (F/RF0)1-U'  as  u>  decreases.  The  effect  of 

wall  temperature  Tw  is  found  to  be  small.  At  lower 

where  x  is  the  distance  of  the  point  considered  “a11  temperatures  the  flow  reversal  is  delayed  to 

from  the  initiation  position.  Figure  45  shows  the  some  extent,  compared  to  the  higher  wall  tempera- 

variation  of  heat  transfer  to  the  wall  with  time  t,  ture. 
at  different  positions  for  plane  flow.  Also  shown 

in  Fig.  45  are  the  experimental  results  from  Ref.  To  test  the  analysis,  the  results  of  heat 

28.  The  agreement  between  the  analytical-numerical  transfer  to  the  wall  behind  a  planar  detonation 

results  and  the  experimental  data  is  satisfactory.  wave  ®°ving  into  a  stationary  mixture  of  2Hj*0, 

A  better  relation  for  the  temperature-viscosity  were  computed  and  the  variation  of  heat  transfer 

might  even  improve  the  agreement.  with  tine  was  compared  with  experimental  data  at 

different  positions  from  the  initiation  point. 

To  check  the  influence  of  the  initial  conditions  Very  good  agreement  was  obtained  between  the 
on  heat  transfer,  calculations  were  also  done  for  present  numerical  results  and  the  experimental 

the  variation  of  the  average  heat  transfer  over  data.  This  lends  confidence  to  the  present  ana- 

given  periods  of  time  for  different  pressures  p  ,  lytical-numerical  work  which  is  less  restrictive 

and  is  shown  in  Fig.  46.  The  crosses,  closed  and  in  it*  assumptions  and  models  the  actual  flow  more 

open  circles  shown  in  Fig.  46  represent  similar  realistically  than  has  been  done  previously, 

experimental  results  from  Ref.  29.  It  can  be  seen 

from  Fig.  46  that  the  relation  between  heat  transfer  The  boundary  layer  behind  a  detonation  wave 

and  initial  pressure  is  not  linear  and  does  not  may  undergo  transition  to  a  turbulent  one.  Conse- 

justify  the  conclusion  in  Ref.  29  that  it  is  linear.  quently,  a  complete  boundary-layer  solution  through 

the  transition  and  turbulent  flow  is  very  difficult 
As  shown  in  Ref.  8,  the  structure-independent  and  must  be  done  in  stages,  as  noted  in  the  Intro- 

detonation  inviscid  flow  properties  in  a  2Hj*0j  duction.  For  the  strong  blast  wave,  the  density 

mixture  are  very  insensitive  to  the  initial  condi-  behind  the  wave  front  decreases  so  rapidly  that  the 

tions.  For  example,  the  detonation  speed  changes  inviscid  flow  becomes  very  rarefied  and  the  boundary- 

from  2.59  km/sec  to  3.08  km/sec,  the  temperature  layer  thickness  increases  very  quickly  in  a  short 

from  2934  K  to  4645  K,  and  the  pressure  ratio  from  distance  (say,  0  <  E.  <  0.3)  away  from  the  wave  front. 

16  to  21.8  even  when  the  initial  pressure  changes  Consequently,  the  boundary-layer  concept  may  no 
from  0.01  atm  to  100  atm  at  an  initial  temperature  longer  be  valid  in  the  region  beyond  £  >  0.3.  To 

of  298.15  K.  To  simplify  the  procedure  of  calcula-  illustrate  this  poi.n  several  boundary- layer  profiles 

tion  an  average  value  of  B)  was  taken  in  calculating  have  been  plotted  (Fig.  44)  and  discussed.  The 
the  heat  transfer  with  reasonable  approximation.  present  report  represents  the  first  of  several 

Figure  47  shows  a  comparison  of  the  present  results  stages  as  noted  above  in  the  solution  of  the  boun- 

with  experimental  and  other  analytical  results  from  dary  layer  in  the  UTIAS  implosion  chamber.  Never- 

Ref .  15.  It  can  be  seen  that  present  analysis  is  theless,  the  results  are  applicable  to  many  other 

in  very  good  agreement  with  the  experimental  data  problems. 
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Table  1 


inviscid  Flow  Solution  Behind  a  Strong  Spherical  Blast  Wave 


0=2,  m  =  2/5 

e 

9> 

F 

R 

0.00 

0.8333 

0.8333 

6.0000 

i  0.10 

0.68486 

0.4231 

1.2318 

i  0.20 

0.58374 

0.33714 

0.39251 

0.30 

0.50309 

0.31324 

0.13251 

-1 

0.40 

0.42927 

0.30656 

0.40661x10 

j  o.so 

0.35738 

0.30495 

0. 10284X10*1 
-2 

•  0.60 

0.28596 

0 . 30464 

0.19214x10 

0.70 

0.21472 

0 . 30460 

0.220S5xl0"3 

-4 

0.80 

0.14384 

0.30460 

0.10238x10 

-7 

0.90 

0.075392 

0.30460 

0.40588x10 

Table  2 


Flow  Derivatives  Behind  a  Strong  Spherical  Blast  Wave 


0=2,  m  =  2/5 

« 

Fe 

.  . 

0£^ 

0.00 

1.8056 

13.028 

125.000 

0.10 

1.1956 

2.2153 

15.560 

0.20 

0.87500 

0.58514 

4 . 2368 

0.30 

0.75925 

0.16773 

0 .43694x10* 

_-2 

1.4782 

0.40 

0.7243S 

0.51333 

0.50 

0.71549 

0.9187x10 

0.15468 
.  _-l 

0.60 

0.71326 

0.13726x10  2 

0.36112x10 

0.70 

0.71136 

0.11818xl0'3 

-S 

O.SS426xIO*2 

0.80 

0.70435 

0.36594x10 

.  -8 

0.39075x10  0 

0.90 

0.634  78 

, _ 

0.73131x10 

0.35366x10 

The  derivatives  are  all  negative. 
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Table  3 


Solution  for  lnviscid  Flow  Field  Behind  C-.)  Detonation  Wave 


o  = 

0,  m  = 

1 

0  = 

1 ,  m  = 

1 

O  S 

2 ,  m  = 

1 

c 

F 

R 

«p 

F 

K 

«p 

F 

R 

0.00 

0.4544 

0.47752 

1 .8328 

0.4544 

0.47752 

1.8328 

0.45440 

0.4775 

1 .8328 

0.05 

0.40772 

0.43301 

1.6821 

0.34050 

0.36820 

1.4595 

0.30781 

0.33779 

1.3534 

0. 10 

0.36105 

0.39234 

1.5430 

0.28650 

0.32261 

1.3000 

0.24717 

0.28807 

1 . 1773 

0.15 

0.31437 

0.35527 

1.4146 

0.24088 

0.28807 

' .1773 

0.19961 

0.25350 

1 .0527 

0.20 

0.26769 

0.32151 

1.2962 

0.19917 

0.25966 

1 .0750 

0.15871 

0.22697 

0.9556 

0.25 

0.22101 

0.29077 

1.1870 

0.15984 

0.23552 

0.98700 

0.12222 

0.20576 

0.8769 

0.30 

0.17434 

0.26281 

1.0865 

0.12234 

0.21474 

0.1)1035 

0.089251 

0.18851 

0.8122 

0.3S 

0.12766 

0.23738 

0.99388 

0.086556 

0.19684 

0.-4356 

0.059572 

0.17450 

0.7591 

0.40 

0.08098 

0.21428 

0.90866 

0.052765 

0. 18155 

0.78591 

0.033490 

0.16146 

0.7092 

0.45 

t 

0.034304  j 

0.19329 

0.83026 

0.021685  j 

0.16881 

0.73742 

0.012032 

0.15504 

0.6845 

Table  4 


Derivatives  of  lnviscid  Parameters  Behind  C-J  Detonation  Wave 


5  = 

0,  ra  = 

1 

0  = 

1 ,  m  = 

1 

0  = 

2,  m  = 

1 

f. 

*5 

Fc 

Rc 

Fc 

Rc 

Fc 

Rc 

0.93356 

0.93355 

3.1361 

-00 

153.840) 

— cc 

(53.827) 

-CO 

(180.92) 

—CO 

(75.720) 

(75.702) 

(254.51) 

0.05 

0 .93356 

0.85156 

2.8958 

1.2399 

1.1030 

3.8274 

1.4515 

1.2615 

4.4247 

0.10 

0.93356 

0.77633 

2.6727 

0.97132 

0.77469 

2.7328 

1.0467 

0.8045 

2.8783 

0.15 

0.33356 

0.70735 

2.4655 

0.86560 

0.62075 

2.2208 

0.87356 

0.59808 

2.1741 

0.20 

0.93356 

0.64412 

2.2732 

0.80745 

0.52153 

I .8902 

0.76941 

0.47150 

1.7377 

0.25 

0.93356 

0.58620 

2.0949 

0.76743 

0.44701 

1 .6399 

0.69293 

0.38147 

1 .4232 

0.30 

0.93356 

0.53317 

1.9295 

0.73318 

0.38556 

1 .4309 

0.62663 

0.31085 

1.1725 

0.35 

0.93356 

0.48465 

1.7763 

0.69719 

0.33137 

1.2432 

0.55954 

0.25079 

0.95509 

0.40 

0.93356 

0.44028 

1.6344 

0.65220 

0.28050 

1 . 0630 

0.48058 

0.19507 

0.74903 

0.45 

0.93356 

0.39972 

1.5030 

0.58542 

0.22807 

0.87216 

0.36919  j 

0.13595 

0.52541 

The  derivatives  are  all  negative. 
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Table  S 


Comparison  of  Wall  Derivatives  Obtained  by  Three  Different  Numerical  Methods 


- , 

Vf-°> 

Gn(C.0) 

c 

Series 

Expansion 

Method 

Liu 

and 

Mirels 

Solution 

! 

Present 

Series 

Expansion 

Method 

Liu 

and 

Mirels 

Solution 

1 

Present 

a 

0.66141 

0.66198 

1 

0.66147 

0.89693 

0.89864 

1 

0.89735 

Table  6 

Comparison  of  Wall  Derivatives  of  Liu  and  Mirels  (1980)  and  Present  Work 
for  a  Boundary  Layer  Behind  a  Strong  Blast  Wave 


gn(C.°) 

e 

Liu  and 
Mirels 

Present 

Liu  and 
Mirels 

Present 

0.00 

0.66198 

0.66147 

0.89864 

0.89735 

0.10 

0.99598 

0.99520 

0.71326 

0.70950 

0.20 

1.4215 

1.4210 

0.6506 

0.6486 

0.30 

2.0170 

2.0170 

0.67915 

0.6789 

0.40 

2.8813 

2.8800 

0 . 8664 

0.8692 

0.50 

4.03096 

4.0240 

1.2942 

1 . 2980 

0.60 

5.52177 

5.5090 

1.8290 

1.8400 

0.70 

7.94119 

7.9040 

2.61496 

2.6490 

0.80 

12.7369 

12.570 

4.o957 

4.3220 

0.90 

27 . 2309 

25.26 

9.0235 

9.770 
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Table  7 

Comparison  of  Wall  Derivatives  for  Different  Cases  Behind  a  Chapman-Jouguec  Wave 
iTw  f_  AQP_  Ajl  -  »  0.75,  Pr  -  2.26) 


=  0 

ii  =  0 

a  *  1 

(7  =  0 

O  =  1 

a  =  1 

I  o  =  2 

c  *  1 

f. 

f  (f.,0) 

nn 

8n(t.°) 

f  (*..0) 

on 

,  8n(f-0) 

1  f  (f.,0) 

nn 

8n(f-°)  | 

0.0 

0.4023 

0.6354 

0.4023 

0.6354 

0.4023 

0.6354 

0.4023 

0.6354 

0.05 

0.4149 

0.6506 

0.4421 

0.6883 

0.4673 

0.7269 

0.4800 

0.7482 

0. 10 

0.4292 

0.6699 

0.4586 

0.7191 

0.5145 

0.8043 

0.5275 

0.8369 

0.20 

0.4467 

0.7101 

0.4656 

0.7742 

0.5943 

0.9806 

0.5882 

1.0330 

0.30 

0.4293 

0.  7587 

0.4083 

0.8314 

0.6069 

1.2120 

0.  5190 

1.2830 

0.40 

0.2145 

0.8194 

0.05685 

0.8947 

0.1353 

1 . 53*0 

-0.3740 

1.6190 

0.45 

-0.5084 

0.8537 

-0.9738 

0.9248 

-1.7140 

1.7410 

-3.6520 

1.8320 

0.50 

-0. 1279 

0.8895 

-0.4341 

0.9596 

-0.8039 

1.9860 

-2. 3680 

2.0820 

0.60 

0.2235 

0.9310 

0.04352 

0.9994 

0.2114 

2.6110 

-0.9636 

2.7320 

0.70 

0.4271 

0.9585 

0.3287 

1.0270 

1.2560 

3.6320 

0. 3380 

3. 7980 

0.80 

0.5605 

0.9783 

0. 5208 

1.04  70 

2.9270 

5.  7040 

2.2290 

5.9640 

0.90 

0.6513 

0.9933 

0.6577 

1.0620 

7.6860 

12.1500 

7.3090 

12.690 

0.95  | 

0.6836 

0.9996 

0.7111 

1.0690 

j  17.360 

|  24.350 

1  17.430 

.  25.250 

Table  8 


Ratios  of  Boundary  Layer  Thicknesses  for  Different  Prandtl  Number,  Pr 


f. 

Pr=0.72 

Pr=l .0 

Pr=l . 5 

Pr-2.256 

V5T 

V#T 

(5  /fi 

V  T 

V7  T 

0.00 

0.10 

0.8644 

1.000 

1.377 

1.791 

0.20 

0.7695 

1.000 

1.303 

1.749 

0.30 

0.6486 

1.000 

1.364 

1.702 

0.40 

0.  3456 

1.454 

1.833 

0.50 

0.2129 

1.000 

1.541 

1.862 

0.60 

0.2676 

1.490 

1.858 

0.70 

0.2677 

1.522 

1.812 

0.80 

0.4186 

1.514 

1.941 

0.90 

0.6416 

1.000 

1.278 

2.263 

0.95 

0.7407 

1.000 

1.5120 

2.158 
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Table  9 


Wall  lVrivati  ;  for  Different  Flow  Conditions  behind  a  Chaproan-Jouguet  Wave 

<n  =  2,  a  =  1) 


! 

T  -300  K 

w 

'  T  =300  K 

'  T  =300  K 

w 

! 

T  = 
w 

■0.0  J 

o=0.  75 

I  Pr=2 . 256 

1 

1  u>=0. 75 

_  1 _ 

Pr=0. 72 

u>=1.00 

Pr=2. 256 

u=l .  0 

!  Pr=2. 256  ■ 

"0> 

!  £,n(;’,0) 

gn(C,0) 

f  (5,0) 
no 

gn(C.0) 

f  (5,0) 

nn 

1  erU>°> 

0.0 

0.4023  ! 

0.6354 

0.4171 

0. 3352 

0.7281 

1.115 

0.7281 

1  1.195 

0.05 

0.4800 

0.7482 

j  0.5296 

0.4038 

0.8675 

1.309 

0.8774 

1.412 

0. 10 

0.5275 

0.8369 

;  0.6125 

0.4558 

0.9505 

'  1.456 

0.9716 

1.577 

0.20 

0.5882 

1.0330 

!  0.7985 

0.5699 

1.061 

1.785 

1.123 

1.941 

0.  30 

0. 5190 

1.2830 

1.058 

0.7147 

0.9547 

2.199 

1.123 

2.400 

0.40 

-0.3740 

1.6190 

1.567 

0.9091 

-0.5231 

2.751 

0.0631 

3.007 

0.45 

-3.652 

1.8320 

2.540 

1.032 

-0.6002 

3.1010 

i 

-0.4284 

3.  389 

0.  50 

-2.368 

2.0820 

3.138 

1.172 

-4.079 

3.500 

-3.203 

3.821 

0.60  : 

-0.9636 

2.7320 

3.541 

1.525 

-1.834 

4.526 

-1.318 

4.935 

,  0.70  ■ 

0.3380 

3.7980 

4.152 

2.101 

0.2729 

6.161 

0.6444 

6.714 

0.80  : 

2 . 2290 

5.9640 

5.533 

3.67 

3.247 

9.  322 

3.538 

10.161 

0.90 

7.3090 

12.690 

10.15 

6.928 

10.52 

18.08 

10.76 

19.  70 

0.95  j 

17.430 

25.250 

|  19.72 

14.57 

i  22.85 

31.75 

i 

23.04 

34.59 
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FIG.  2  NORMALIZED  INVISC3D  FLOW  PROFILE  .  A  STROx.u  SPHERICAL 

BLAST  WAVE  FOR  DENSITY  R,  PRESSURE  f  AND  FLOW  VELOCITY  <p 
WITH  DISTANCE  | . 


FIG.  3  NORMALIZED  PRESSURE  DERIVATIVE  F|,  VELOCITY  DERIVATIVE  <p . 

AND  DENSITY  DERIVATIVE  R|  FOR  INVISCID  FLOW  BEHIND  A  STRONG 
SPHERICAL  BLAST  WAVE  WITH  DISTANCE  ?  AT  |  =  0,  <?>t  =  1.8056, 
F|  =  13.028,  Rj  =  1.25. 


FIG. 


e 


NORMALIZED  INVISCID  FLOW  VELOCITY  PROFILES  <p  BEHIND  C-J  DETONATION 
WAVES  WITH  DISTANCE  |  FOR  PLANAR  (<f  =  0),  CYLINDRICAL  (<f  =  l)  AND 
HEMISPHERICAL  (cf  =  2)  FLOWS. 


1 
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FIG.  7  NORMALIZED  INVISC1D  FLOW  DERIVATIVES  <P|  FOR  VELOCITY  BEHIND 
C-J  DETONATION  WAVES  WITH  DISTANCE  |  FOR  PLANAR  (a  =  0), 
CYLINDRICAL  (&  »  l)  AND  SPHERICAL  (a  =  2)  FLOWS. 


Fe 


FIG.  9  NORMALIZED  INVISCID  FLOW  DERIVATIVES  Ft  FOR  PRESSURE  BEHIND  C-J 
DETONATION  WAVES  WITH  DISTANCE  g  FOR  PLANAR  (cr  =  0),  CYLINDRICAL 
(a  =  1)  AND  SPHERICAL  (c  =  2)  FLOWS. 


FIG.  12  NORMALIZED  VELOCITY  fT,  AND  TEMPERATURE  g  PROFILES  WITH  DISTANCES 


4.00 


FIG.  13 


NORMALIZED  VELOCITY  PROFILE  IV 
T)  AND  ?  BEHIND  PLANAR  C-J  DETO 


Tw  =  300  K,  a  =  0,  cr  =  0. 


OF  BOUNDARY  LAYER  WITH  DISTANCES 
NATION  WAVE.  Pr  =  2.26,  w  =  0.75, 


5.00 


4.00 


FIG.  15  NORMALIZED  VELOCITY  PROFILES  f„  OF  BOUNDARY  LAYER  WITH  DISTANCES 
r\  AND  |  BEHIND  nYLINDRICAL  C-J  DETONATION  WAVE.  Pr  =  2.26 
w  =  0»75j  Tw  =  300  K,  cr  =  1,  a  =  0. 


*  i—  - 


- 


5.00 


4.00 
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FIG.  17  NORMALIZED  VELOCITY  PROFILES  fV  OF  BOUNDARY  LAYER  WITH  DISTANCES 
T]  AND  |  BEHIND  CYLINDRICAL  C-J  DETONATION  WAVE.  Pr  =  2.26, 
w  =  0.75}  Tw  =  300  K,  a  =  1,  a  =  1. 
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FIG.  21 


VARIATION  OF  NORMALIZED  VELOCITY 
DERIVATIVE  g^i, 0)  WITH  DISTANCE 
Pr  =  2.26,  w  =  0.75,  Tw  =  300  K, 


DERIVATIVE 
|  FOR  PLANAR  C 


(1,0)  AND  TEMPERATURE 
-J  DETONATION  WAVE. 


cr  =  0,  o'  =  0. 


J 


8 


97?(£,0) 


FIG.  22  VARIATION  OF  NORMALIZED  VELOCITY  DERIVATIVES  f«T,(|,0)  AND  TEMPERATURE 
DERIVATIVE  grj(l,0)  WITH  DISTANCE  |  FOR  CYLINDRICAL  C-J  DETONATION 
WAVE.  Pr  =  2.26,  go  =  0.75,  Tw  =  300  K,  <7  =  1,  cr  =  0. 
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FIG.  23  VARIATION  OF  NORMALIZED  VELOCITY  DERIVATIVES  fnT1(|,0)  AND  TEMPERATURE 
DERIVATIVES  gT)(|,0)  WITH  DISTANCE  |  FOR  CYLIND®AL  C-J  DONATION 
WAVE.  Pr  =  2.26,  w  =  0.75,  Tw  -  300  K,  5  =  1,  <r  =  i. 


2.0 


FIG.  2k  VARIATION  OF  NORMALIZED  VELOCITY  DERIVATIVES  fnn(l,0)  AND  TEMPERATURE 
DERIVATIVES  g-UjO)  WITH  DISTANCE  |  FOR  SPHERICAL  C-.T  DETONATION 
WAVE.  Pr  =  2;26,  w  =  0.75,  Tw  =  300  K,  a  =  2,  a  =  1. 


1L. 


5.00 


FIG.  25 


NORMALIZED  VELOCITY  PROFILES 
I  AND  T)  BEHIND  SPHERICAL  C-J 
Tw  »  300  K,  5  =  2,  <r  =  I. 


fn  OF  BOUNDARY  LAYER  WITH  DISTANCES 
DETONATION  WAVE.  Pr  =  1.5,  w  =  0.7 


PROFILE  f3  OF  BOUNDARY  LAYER  WITH  DISTANCES 
ICAL  C-J  DETONATION  WAVE.  Pr  =  1,  w  =  0.75, 


FIG.  2d  NORMALIZED  TEMPERATURE  PROFILER  t'  OF  BOUNDARY  LAYER  WITH 
DISTANCES  tj  AND  £  BEHIND  SPHERICAL  C-.T  DETONATION  WAVE. 
Pr  =  1,  w  -  0 . 7‘> ,  Tw  -  300  K,  cr  =  2,  a  =  1. 
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FIG.  30  NORMALIZED  TEMPERATURE  PROFILES  g  OF  BOUNDARY  LAYER 
T]  AND  |  BEHIND  SPHERICAL  C-J  DETONATION  WAVE  Pr  = 
Tw  =  300  K,  ct  =  2,  ct  =  1. 


WITH  DISTANCES 

0.72,  u>  =  0.75 


FIG.  31  EFFECT  OF  Pr  ON  NORMALIZED  VELOCITY  DERIVATIVES  fnnU,0)  WITH 

DISTANCE  |  FOR  SPHERICAL  C-J  DETONATION  WAVE,  u>  =  0.75,  Tw  =  300  K 
<7  =  2,  cr  =  L.  5 


FIG.  35  NORMALIZED  TEMPERATURE  PROFILES  g  OF  BOUNDARY  LAYER  WITH 
DISTANCES  t,  AND  |  BEHIND  C-J  SPHERICAL  DETONATION  WAVE 
Pr  =  2.25,  u  =  L,  Tw  =  300  K,  a  =  2,  a  =  1. 


5.00 


4.00 


FIG.  37  NORMALIZED  TEMPERATURE  PROFILES  g  OF  BOUNDARY  LAYER  WITH 
DISTANCES  T)  AND  |  BEHIND  SPHERICAL  C-J  DETONATION  WAVE. 
Pr  =  2.26,  U)  =  1,  Tw  =  0  K,  a  =  2,  a  =  1. 
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FIG.  38  VARIATION  OF  NORMALIZED  VELOCITY  DERIVATIVES  f,r(Lo)  WITH  DISTANCF 
5  FOR  SPHERICAL  C-J  DETONATION  WAVE.  Pr  =  2.2^,‘  Tw  =  300  K  a  ---  2 
a  =  I  (for  DIFFERENT  VISCOUS  EXPONENT,  w).  ’ 


(£  ,o) 


r  — ^ 


FIG.  39 


VARIATION  OF  NORMALIZED  VELOCITY  DERIVATIVES  fn, 
I  FOR  SPHERICAL  C-J  DETONATION  WAVE.  Pr  =  2.26, 
a  =  I  (FOR  DIFFERENT  WALL  TEMPERATURE,  Tw). 


U,0)  WITH  DISTANCE 
«  =  0.79,  a  =  2, 


FIG.  *+0  VARIATION  OF  NORMALIZED  TEMPERATURE  DERIVATIVES  g  (t  0)  WITH 
DISTANCE  ?  FOR  SPHERICAL  C-J  DETONATION  WAVE.  Pr  =  2.26, 

Tw  =  300  K,  a  =  2,  a  =  1  (FOR  DIFFERENT  VISCOUS  EXPONENT  1  c>) 


(£» 


FIG.  Ul  VARIATION  OF  NORMALIZED  TEMPERATURE  DERIVATIVE  g„(  f,  o)  WITH 
DISTANCE  |_FOR  SPHERICAL  C-J  DETONATION  WAVE.  Pr  =*2.26, 
w  =  0.75,  5  -  2,  a  =  1  (FOR  DIFFERENT  WALL  TEMPERATURE,  Tw) . 
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FIG.  42  VARIATION  OF  THERMAL  BOUNDARY -LAYER  THICKNESS  &r  WITH  DISTANC1 
%  FOR  UTIAS  HEMISPHERICAL  IMPLOSION  CHAMBER.  Pr  =  2.2b  u>  = 
Tw  -  300  K,  Xs  =  10  cm,  Us  =  2981  m/sec,  a  =  2,  a  =  1. 
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FIG.  43  VARIATION  OF  THERMAL  BOUNDARY -LAY HR  THICKNESS  4)T  WIT 
AT  ORIGIN  FOR  UTIAS  HEMISPHERICAL  IMPLOSION  CHAMBER. 
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FIG.  U6  COMPARISON  OF  EXPERIMENTAL  (  Q,  •  AND  o  FROM  REF.  29)  AND 

ANALYTICAL  ( - , - AND - ,  PRESENT  RESULTS)  VARIATIONS 

OF  AVERAGED  HEAT  TRANSFER,  qw  (BTU/ft2  sec)  TO  THE  WALL, 

OVER  50,  100  AND  150  ns,  RESPECTIVELY,  WITH  INITIAL  PRESSURE 
Poo  (TQRR)  BEHIND  A  HYDROGEN- OXYGEN  GASEOUS  DETONATION  WAVE. 
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APPENDIX  A 


DERIVATION  OF  TRANSFORMED  BOUNDARY -LAYER  EQUATIONS 


The  unsteady,  compressible  laminar  boundary- layer 
equations  for  a  perfect  gas  are: 


31  mt 

3x  =  '  ~u  ‘ 

s 


(AM) 


Continuity: 

3p  .  1  3(pux°)  .  3(pv) 

3t  xo  57~  3y  " 


(Al) 


x°  i 
r- 


.At2mo*2u..(ni-l)»l[|l/2 


IA1S) 


Momentum : 


P 


Du 

Dt 


K  +  ± 

3x  +  3y 


Energy : 


P 


Dh 

Dt 


°Pe 

Dt 


(  v  ah 

3u 

i  Pr  3y  J 

.  3y  . 

(A2) 


(A3) 


an 

3t 


3f, 

at 


( 1  -f, )  mt 


(Alt>) 


State : 

P  =  [(-y-D/vJph  (A4) 

To  eliminate  the  explicit  dependence  on  time, 

Mirels  and  llamman  [Ref.  18]  devised  a  similarity 
transformation  for  the  lamina  boundary-layer  induced 
by  a  strong  blast  wave.  Only  a  small  modification 
to  their  transformation  was  needed  to  make  it  applic¬ 
able  either  to  blast  waves  or  for  detonation  waves. 
The  modified  transformed  coordinates  become 


Introducing  a  scalar  stream  function  defined  by 


px 


30- 
o  3y 


(A18) 


px 


3^  ,  3 
37  3t 


td> 


(A19) 


T  =  t 


(AS) 


jAt2ma+2oi(m-l)  +  l^j  1/2 


(A6) 


(A7) 


For  the  new  variables  £,  n  and  T,  the  derivatives 
with  respect  to  x,  y  and  t  become 


j_  3r  _3_  3£j_  3n  _3_ 
at  at  3t  +  at  3f  *  at  3n 

a  -  as  _a_  an  _a_ 

3x  3x  at;  *  17  3n 


(A8) 

(A9) 


a  an  a 
W  =  5y  3n 


(A10) 


substituting  Eqs.  (A13)  to  (A19)  into  Eqs.  (A8)  to 
(A10) ,  and  adding  Eqs.  (A10)  to  (A9) ,  the  following 
derivatives  are  then  obtained 


2mo->2u)(m-l )  +  lf  j 


_ a_ 

172  an 


(l-f)mt 


1  3 

31, 


2o+a(2« -1 ) 


(A20) 


_a_  _a_  f!!_3ndcfi^i_  aj;  _a_  ] 

3x  *  3y  pxo  3y  dx  (  3n  3C  3 C  Bn  | 


jAt2mo+2w(m-l)*l^j 1/2  8n 
Introduce  the  similarity  parameters 


(A21 ) 


For  the  inviscid  flow,  we  are  interested  in  self¬ 
similar  motions  in  the  form  of 


Then 

Using  Eqs. 
follows: 


*s  -  ct 


(All) 


u  -  cmt 
s 


i-l 


(A12) 

(All)  and  (A12)  the  derivatives  are  as 


3t 

at 


l 


F(E)  * 


««  -  if 

s 


(A22) 


( A23) 


(A13) 


(A24) 


and  the  temperature  ratio 


and  note  that 


Te  °»  2  F 

T  =  b  p  Us  R 


(A25) 


2o 


where 


b  = 


Y„(Y-1)C 


for  the  inviscid  flow.  The  dimensionless  stream 
function  is  defined  as 


f(E,n)  -  */{ue[AT2mo+2“(,n-1)+1  ei1/2} 
The  total  derivative  is  defined  by 


(A26) 


D  3  3  3 

P  Dt  -  0  +  PU  3i  *  pV  37 


(A27) 


P  ^  =  -pmt 


•1  n 

'  L  2 


2E 


[2o+a(2w-l)  -  (<pfc+<P£f) 

)  _L 

T)J  3E 


+  *  ^f>J  i  -  -k  } *  p  #r  <A28> 


For  the  momentum  equation,  Eq.  (A2) ,  the  transform¬ 
ation  can  be  done  term  by  term.  For  p(Du/Dt),  we 
have 


Du 

p  Dt  =  “pmt 


v{[? 


(2o*a(2oj-l)  -  (?f?+«P^f) 


2T  <"-«>]  fnn  -  5  )  fn  '  afn} 


(A29) 


For  3pe/3x,  using  Eqs.  (A22)  to  (A24),  the  following 
expression  is  obtained: 


3pe  -1  g  FC 

=  ^s^  R  ^ 


(A30) 


where  g  =  h/he  *  T/Te  is  the  enthalpy  or  temperature 
ratio.  A  new  dependent  variable  is  defined  as 
follows: 


B 


pp 


With  the  transformed  coordinates,  E,  n,  the  term 
<*/3y[u(3u/3y)]  can  be  expressed  in  the  following 
form 


2o  p  „ 

X  u  p  ii  9 

Jr  (p  ly  )  *  [At*^2(m-'lW*ic]  (Bfnn)n  (A32) 


Define 

2C2(0+u))  -2-1  w.  r  cm  T 

A  *  p.  Lb  p„  F°J 


jAt2raa*2(m-l)iDflf  jt-l 


d-n20  C2(0+w)  -> 

~r 


Au 


2u) 


Substituting  Eqs.  (A29) ,  (A30)  and  (A32)  into  the 
momentum  equation  and  dividing  both  sides  by 
pmusqpt~  ,  the  following  equations  are  obtained: 


(l-E)20  R^fF/F  )W(Bf  )  ♦  (n-Vf)f 
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f  .  fne ) 
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(A33) 


for  the  momentum  equation. 


Substituting  Eqs.  (A22)  to  (A24)  and  (A26J  into  (A20) 
and  (A21),  then  Eq.  (A27)  becomes 


The  same  procedure  can  be  applied  to  the 
energy  equation,  Eq.  (A3),  term  by  term: 
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Noting  that 


e  e 
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=  R^Lb^ussJ  ,A38) 
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the  transformed  energy  equation  becomes 
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APPENDIX  B 


THERMAL  CONDUCTIVITY  AND  VISCOSITY  FOR  MULTICOMPONENT  GAS  MIXTURES 


The  viscosity  u  and  thermal  conductivity  K 
contained  in  the  momentum  and  energy  equations  are 
estimated  by  using  the  following  semi-empirical 
relations  given  in  Ref.  30. 


n  N.y. 
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V  =  I 


i=l 


y  n.0.. 

j=i  >  12 


for  viscosity  of  gas  mixture,  and 

n  N.K. 

K  -  l  —■ 1  1 

i=l  " 


l  N.d 
j=l  J 


(Bl) 


(B2) 


for  thermal  conductivity  in  which 
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and  n  is  the  number  of  chemical  species  in  the 
mixture;  Nj  and  Nj  are  the  mole  fractions  of  species 
i  and  j;  Uj  and  U;  are  the  viscosities  of  species 
i  and  j  at  the  system  temperature  and  pressure;  and 
Mj  and  Mi  are  the  corresponding  molecular  weights. 
Note  that  d>jj  is  dimensionless  and  when  i  =  j , 

=  1.  Basically,  Eqs.  (Bl)  and  (B2)  are  adequate 
only  at  low  density.  As  shown  in  Figs.  Bl  and  B2, 
the  viscosity  and  thermal  conductivity  of  a  gas 
approach  a  definite  limit  (the  low-density  limit)  as 
the  temperature  reaches  a  very  high  value  at  a  given 
pressure.  The  temperature  of  the  gas  mixture  behind 
a  C-J  wave  is  about  4000  K.  Taking  water-steam  as 
an  example,  the  critical  temperature  is  400  K  and 
the  critical  pressure  is  1  atm.  In  the  region 
behind  a  C-J  wave,  the  average  pressure  is  less  than 
10  atm.  Hence,  its  reduced  pressure  is  less  than 
10  and  its  reduces  temperature  is  10.  From  Figs. 

Bl  and  B2,  it  can  be  concluded  that  the  low-density 
limit  is  still  a  good  approximation. 


APPENDIX  C 


BOUNDARY  LAYER  THICKNESS  CLOSE  TO  THE  ORIGIN 


As  pointed  out  before,  the  origin  is  singular 
for  the  boundary -layer  equations.  Near  the  origin, 
the  solution  becomes  divergent.  Therefore  it  is 
necessary  to  find  an  alternate  way  of  evaluating 
the  thermal  boundary-layer  thickness  at  origin. 
Since  the  gas  velocity  is  very  small  in  the  middle 
of  the  wave  behind  a  C-J  detonation,  the  boundary- 
layer  equations  reduce  to 


3h  3 

3t  * 


u  3h 

Pr  3y 


(Cl) 


l  3T 
a  3t 


(C2) 


where  a  H  K/pC_.  The  corresponding  boundary  condi¬ 
tions  are 


y  =  0,  T  =  T 

w 

y  -*  ®,  T  ■  T 

e 


(C3) 


Assuming  that  the  Pr  and  Cp  are  constant,  Eq.  (Cl)  The  UTIAS  implosion  chamber  is  made  of  steel, 

can  be  rewritten  in  the  following  form  Because  the  thermal  conductivity  of  steel  is  much 


greater  than  that  of  the  gas  and  the  flow  duration 
of  interest  is  very  short,  it  is  reasonable  to 
assume  that  Tw  =  constant,  that  the  gas  thickness 
is  infinite  in  the  y-direction,  its  temperature  is 
uniform  in  the  x-direction,  and  the  natural  convec¬ 
tion  heat  transfer  is  negligible.  Definining  the 
temperature  difference. 


T  =  T  -  T 


For  convention,  we  define  the  thermal  boundary-layer 
thickness  as  =  ye  at  which 

T(ye.t) 


According  to  the  previous  results  for  the  inviscid 
flow,  the  temperature  outside  the  thermal  boundary 
layer  can  be  expressed  as 


Eq.  (C2)  can  be  expressed  as 
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The  boundary  conditions  of  Eq.  (C3)  become 


y  =  0,  T  =  T  =  0 


where  T  =  298.15  K,  b  =  1.2594,  F  =  0.051135, 

R  =  0.669363,  us  =  2981.5  x  102  cm/sec,  then  the 
temperature  is 


Tg  =  12.25  Ta 


Using  these  data,  the  thermal  boundary-layer  thick¬ 
ness  is  obtained  approximately  by 


y  •*  <»,  T  =  T=  T-  T 
e  e  w 


6^.  =  y^  =  1 .8  v^at 


The  solution  of  Eq.  CS)  is  [Ref.  31) 


T  =  (T Jerf 

/4at 


Using  the  relation  of  Eq.  (C4),  the  temperature 
within  the  thermal  boundary  layer  is 

T(y ,t)  =  (T  -T  )erf  ♦  T  (C. 

e  w  /4ot 


According  to  the  definition  of  a  =  K/pCp,  a  simple 
relation  holds: 


<5,j.  *  0.94  /t  cm 


where  t  is  given  in  psec. 


APPENDIX  D 


DERIVATION  OF  HEAT  TRANSFER  EXPRESSION 


According  to  the  definition  of  the  heat  transfer 
rate  to  the  wall,  we  can  express  it  as  follows: 
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Using  the  transformed  coordinate  it  then  becomes 
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Putting  Fqs.  (D3)-(D8)  into  (D2)  and  after  rearrang- 
ine,  we  can  »et 
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If  £,  anti  Re  are  replaced  by 


u  t 
s 


X  ♦  U  t 

s 


Re 


we  can  get 
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Finally  the  heat  transfer  rate  to  the  wall  can  be 
expressed  as 
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FIG.  B1 


VARIATION  OF  NORMALIZED  VISCOSITY  |ir  =  WITH  PRESS! FRF 
AM.  TEMPERATURE  Tr  -  T/Te,  Pr  ,  p/p".  (raOM  R®! 
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FIG.  B2  VARIATION  OF  CONDUCTIVITY  Kr  -  K/K0  WITH  PRESSURE  lr 
AND  TEMPERATURE  Tr  =  T/Te  (FROM  REF.  30). 
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